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ABSTRACT 



From previous studies of the effect of primordial magnetic fields on early 
structure formation, we know that the presence of primordial magnetic fields 
during early structure formation could induce more perturbations at small scales 
(at present 1-10 /i _1 Mpc) as compared to the usual ACDM theory. Matter power 
spectrum over these scales are effectively probed by cosmological observables such 
as shear correlation and Lya clouds, In this paper we discuss the implications 
of primordial magnetic fields on the distribution of Lya clouds. We simulate 
the line of sight density fluctuation including the contribution coming from the 
primordial magnetic fields. We compute the evolution of Lya opacity for this 
case and compare our theoretical estimates of Lya opacity with the existing data 
to constrain the parameters of the primordial magnetic fields. We also discuss 
the case when the two density fields are correlated. Our analysis yields an upper 
bounds of roughly 0.3-0.6 nG on the magnetic field strength for a range of nearly 
scale invariant models, corresponding to magnetic field power spectrum index 
n ~ —3. 

Subject headings: Cosmology: primordial magnetic field, Lya clouds, structure 
formation 
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1. Introduction 



In the past 10 years, cosmological weak lensing and the study of Lya clouds in the 
redshift range 2 < z < 5 have emerged as reliable methods to precisely determine the 
matter power spectrum on scales below 10h _1 Mpc. In particular, these methods can 
estimate the matter power spectrum at small scales which are not directly accessible to 
other methods e.g. galaxy sury e ys (for detai l s anc 



( 120081 ) 



Hoekstra &: Jain 



(120081 ): 



Croft et al. 



( 11998 



urth e r refe rences see e.g. 



Munshi et al. 



1999 



2002|)). 



Large scale magnetic fields been observed in galaxies a nd clusters of g alaxies with the 



Widrowl (120021)). There is also 



19891 ). These 



coherence lengths up to ~ 10-100 kpc (for a review see e.g. 
some evidence of coherent magnetic fields up to super-cluster scales ( jKim et al. 
fields play an important role in various astrophysical processes. Still little is known about 
the origin of large scale cosmic magnetic fields, and their role in the evolutionary history 
of the universe. These fields could have or i ginate d from dynamo amplification of ve r y tiny 



seed magnetic fields ~ 10 20 G (e . 



Parker! ( 119791 ); IZeldovich. Ruzmaikin fc Sokolovl ( 119831 ); 



Ruzmaikin. Sokolov fc Shukurovl (I19881 )). It is also possible that much l arger primordial 



magnetic fie 



1988 



Ratra 



ds (^ 10 9 G) were generated during the inflationary phase (ITurner fc Widrow 



19921 ) and the large scale magnetic field observed today are the relics of these 



fields. These large scale primordial fields could be directly detected by upcoming and future 
radio interferometers such as LOFAR and SKA (for details see e.g. www.lofar.org and 
www. skatelescope . org/pages/page_sciencegen.htm). 

The impact of large scale primordial magnetic field s on CMBR temperature and 



polar i zation anisotropies has been stud i ed in detail ( e.g. 



2002h : 



Seshadri fc Subramanian (2001) 



(2005); 



Tashiro &; Sugivamal ( 20061); 



Sethi &; Subramanian! (12009) 



Giovannini &; Kunzd (120081 ); 



Mack et al 



Su 



(2002); 



jramanian fc Barrow! ( 



Sethi. Hai mam Pandevl (120 10) 



Yamazaki et al 



(120081 ): 



Sethi fc Subramanian! ( 120051); 



Lewis! ( 20041): 



1998B 



Gopa 



Sethi et al 



t Sethi 



(120081): 



Kahniashvili fc Ratra 



Seshadri fc Subramanian 



(|2005l ): 



(120091 ) 



Trivedi et al.l ( 120101 . |2012| )). IWassermanl ( 119781 ) demonstrated that primordial magnetic 
fields can induce density perturbations in the post-recombination universe. Further work 
along these lines have investigated the impact of this effect for the formation of first 
structures, reionization of the universe, and the signal from redshifted HI line from the 



epoch of reionization (e.g 



Kim et al 



1996 



Gopal fc Sethi 



2003 



Sethi fc Subramanian 



2005 



Tashiro &: Sugiyama 



2006 



Schliecher. Banerjee. Klessen 



2009 



Sethi &: Subramanian 



20091 ). The matter power spectrum induced by primordial magnetic fields can dominate 
the matter power spectrum of the standard ACDM model at small scales. Probes such 
as cosmological weak gravitational lensing can directly probe this difference and therefore 
reveal the prese n ce of primordial fields and put additional constraint on their strength 



Pandev fc Sethil tom . 



In this paper we attempt to constrain primordial magnetic fields within the framework 
of the distribution of Lya clouds in the IGM in the redshift range 2 < z < 5. These 
clouds have been show n to originate in the mildly non-linear density regions of the IGM 



( jCen fc Ostrikerl ( 119941 ) ) . This has allowed development of d etailed semi-analytic met 



rods 



to understand the observed p roperties of these clouds (e.g. 



Bi et al. 



1995 



Hui et al 



1997 



Choudhury et al. 



2001a 



b|). Adopting a semi-analytic approach, we simulate density 
fluctuation along the line of sight, including the contribution of matter perturbations 
induced by these magnetic fields. We compute effective Lya opacity of the IGM for this 
computed Lya cloud distrib ution and compare our results with the existing data (e.g. 



Faucher-Giguere et al. 



20081 ). 



In the previous analyses, the density perturbations induced by magnetic fields 
are assumed to b e uncorrelated to t he de nsity field generated by the usual ACDM 



model. Recently, 



Caldwell fc Motta 



2011 



showed that if the conformal invariance of 
electromagnetism is broken during the inflation and thus produced the primordial magnetic 
fields, these magnetic fields may be correlated with the primordial density perturbations. 
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In our analysis we have made a separate case for such fields. 

Throughout the paper, we used flat (k=0) ACDM universe with Q m = 0.24, fib = 
0.044, h = 0.73 and a s = 0.77. 



2. Primordial Magnetic Field And Induced Matter Power Spectrum 



In the primordial theory of the magnetic fields, it is postulated that the large 
scale primordial magnetic fields of strengths ~ 10~ 9 G were present in the very early 
universe; these fields could have originated during the inflation. They are assumed to 
be tangled magnetic fields which can be characterized by a power-law power spectrum: 
M(k) = Ak n . In the prerecombination era, the magnetic fields are dissipated at scales 
below a scale corresponding to fc M ^ 200 x (10~ 9 G/_B e ff) (IJedamzik. Katalinic. fc Olinto 



1998 



Subramanian fc Barrow 



1998AI . e.g.). Here B c g is the RMS at this cut-off scale for a 
given value of the spectral index, n. Another possible normalization, commonly used in the 
literature, is the value of RMS at k = 1 Mpc^ 1 , B$. These two normalizations are related 
as: B e g = Bokmst*^ 2 ■ It is possible to present results using either of the pair {B c g,n} or 
{B Q ,n}. 



The PMF induced matter perturbations grow in the post recombination era by 
gravitational instability. The matter power spectrum of these perturbations is given by 
P(k) oc k 2n+1 , for n < —1.5, the range of spectral indices we consider here fjWasserman 



1978; 



Kim et al. 



1996 



Gopal fc Sethi 



20031). 



Magnetic field induced matter perturbations can only grow fo r scales above the magneti c 



field Jean's length: fcj ~ 15 x (10 9 G/B cff ) (e.g. 



Kim et al. 



1996 



Kahniashvili et al. 



20101). 



The dissipation of tangled magnetic field in the post-recombinati o n era also results in an 



increase in the thermal Jeans length fjSethi fc Subramanian 



2005 



Sethi et al 



20081 ). For 
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most of the range of magnetic field strengths and the physical setting (Lyman-a clouds 
at a temperature of ~ 10 4 K) considered here, the scale corresponding to fcj generally are 
comparable to or smaller than the thermal Jeans length. 

For our computation, we need to know the time evolution of the matter power spectrum 
induced by tangled magnetic fields. It can be shown that the domin ant growing mode in 



this case has the same time dependence as the ACDM model (see e.g. iGopal fc Seth.il (120031 ) 
and references therein) 



3. The Simulation: Density Fluctuation Along The Line of Sight: Distribution 

of Lya Clouds 

We describe a brief outline of the numerical simulation in this section. Hydrodynamical 
simulations have shown that Ly-a clouds are mildly non-linear (5 < 10) regions of the IGM 
at high redshifts. This allows one to analytically derive important observables from the 
Lyman-a clouds semi-analytically, in terms of a few parameters denoting the ionization, 
thermal, and dynamical state of the clouds. 



Here we have closely followed the semi-analytic prescription given in 



Bi &: Davidson 



(|1997l). In this paper we have considered two cases of primordial magnetic field induced 
matter perturbations : (a) pure ACDM matter perturbations and primordial magnetic field 
(PMF) induced matter perturbations are uncorrelated (b) those two are correlated. In both 
cases we compute two separate line of sight density (and velocity) fields each corresponding 
to a single kind of matter perturbations. In the former case, these fields are drawn from 
different realizations and in the latter the fields are generated from the same realization. 
We add these two density (and velocity) fields to get the final density (and velocity) fields 
in the IGM. To simulate line of sight IGM density and velocity fields for a given three 
dimensional matter power spectrum (inflationary/PMF induced), first we calculate the 
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corresponding three dimensional baryon power spectrum, which corresponds to the original 
power spectrum smoothed over the scales below the larger of the thermal or Magnetic Jeans 
scale Xb 

^ {k ' Z) ~ [l+x>(z)&] [ ) 



where 



x b (z) 



27& B T m (z) 



1/2 



H \_3fim p Q, m (l + z) 

then we compute one dimensional baryon (density, velocity and density-velocity) power 
spectra, which will be used in the further computation. We note here that the relevant 
scale of smoothing for the range of magnetic field values and the IGM temperatures we 
consider is thermal Jeans scale and not the magnetic Jeans scale. The one-dimensional 
power spectra can be computed using the following relations 

1 f°° 



(2) 



Pk'(k,z) = — dk'k'Pg>{k,z) (3) 

27T J\ k \ 

P?\k,z)=a\z)e±J^Pg\k,z) (4) 

i POO ri / 

P^{k,z) = ia{z)k- ^rPg\k,z) (5) 

where a is the scale factor. 

The density (5o(k,z)) and velocity (v(k,z)) fields in one dimension are two correlated 
Gaussian random fields (the correlation is given by the density- velocity power spectrum), 
we use inverse Gram-Schmidt procedure to simulate them in terms of two independent 
Gaussian random fields w(k) and u(k) of power spectra respectively P w (k) and P n (k) 



6 (k,z) = D(z)[u(k)+w(k)} (6) 



v(k, z) = F(z)iak/3(k, z)w(k, z) 
7 



(7) 



where D(z) and F(z) are the linear density and velocity growth factors. Functions (3(k), 
P w (k) and P u (k) are function of P^\k), 



Pj 3) /k' 3 dk' 
\k\ P^/k'dk' 



0(M = I * (3 / (8) 



1*1 

1 f°° 

P u (k) = — / Pg\k')k'dk' - P w (k) (10) 

We compute 8o(k,z) and v(k,z) for both kinds of perturbations separately, the 
corresponding 5b(x,z) and v(x, z) is computed by using Fourier transforms. And then we 
add the contribution from both the kinds together (8&(x,z) = <5g fl (x, z) + 8^ (x, z) and 
v(x, z) = v m{l (x, z) + v pmi (x, z)) to get the final combined line of sight density and velocity 
fields. To compute one dimensional density field for the pmf induced perturbations we use 
the three-dimensional matter power spectrum (e.g. Gopal & Sethi, 2003); for inflationary 
perturbations we use the standard ACDM power spectrum. 

For our computations we have generated the density and velocity fields for 25 redshift 
bins between the redshifts to 5. In each bin we have 2 14 points resolving the Jeans scale 
by at least a factor of 4. The cutoff scale (Jeans scale, Xf,) is the larger of the thermal Jeans 
length and the magnetic Jeans length. 

To take into account the non-linearity of d ensity perturba t ions in the IGM we use 



lognormal distribution of the IGM density field 
density of baryons in IGM is taken to be, 



Bi fc Davidson! (119971 ). thus the number 



n B {x,z) = Ae 5B{x ' z) (11) 



S 



A is a constant which can be determined using following relation, 

(n B (x,z))=n (z) = A(e 5 ^) (12) 

since Sb{x,z) is a Gaussian random variable, 

( e M^)} = e <«l(*,*)> ( 13 ) 

thus 

n B (x, z) = no (z)e^ z) - {s ^ z))) (14) 
where no (z) is the background baryon number density given by, 

n,{z) = ^{l + zf (15) 

4. Calculation of Lya opacity 

The optical depth r is given by 

J n m (t)a a (£) dt (16) 



r iv) 



where nni is number density of neutral hydrogen in the IGM, v is the observed frequency, 
which is related to redshift z by z = [y a l v ) — 1, is the Lya frequency at rest. The 
absorption cross section a a is given by, 

o"a = 7-^^ a, —7 1- T ( 17 ) 



fe^Tr V ' bv a b 

where parameter b = (2kT/m p ) 1 ^ 2 is the velocity dispersion and v(x) is the peculiar velocity 
field, a = 27re 2 v a /3m e c 3 b = 4.8548 x 10~ 8 /6, I a = 4.45 x 10" 18 cm- 2 and V is the Voigt 
function. 

The number density of neutral hydrogen, Uhi can be computed by solving ionization 
equilibrium equation, 

9 



nm(x,z) 



at[T(x, z)]n B (x, z) 



18) 



a[T(x,z)] + T ci [T(x,z)} + J(z)/[n e n B (x, z)] 
where T(x,z) is given by T(x,z) = T (z)[n B (x, z)/n (z)} 1 ^ 1 where T (z) is the temperature 
of the IGM at the mean density and 7 is the polytropic index for the IGM; 7 captures 
the dynamical state of the IGM gas which gives rise to the observed Lyman-a absorption. 
These parameters are likely to take values in the ranges 4000 < To < 15,000 K and 1.3 
^ 7 ^ 1-6 (Hui & Gnedin, 1997). a(T), T ci (T) and J(z) are recombination rate, collisional 
ionization rate and photo ionization rates in the IGM. For temperature T ~ lO^K, the 
combination of these effects yields (Croft et al. 1998), 



t(z) oc n D T 



-0.7 



A(n B /n Q ) 



2-0.7(7-1) 



A = 0.946 



l + z 



4 J Vo.0125; \io 4 k) Vio 12 ^- 1 



-0.7 



J 



-1 r 



H{z) 

Hn 



-1 



(19) 



To compare with the data we have computed effective optical depth which is the 
observable quantity in the form of decrease in observed flux (F oc e~ T ) and is given by, 



Teg{z) = - log [(exp(-r))] 



(20) 



The data which we have used for comparison with simulation results has been obtained 



using high resolution spectral observations like HIRE S , ES 



in the range of 6-44 km sec 1 (IFaucher-Giguere et al. 



( 12008 



and MIKE having FWHM 
)), which resolve the Jeans 



scales over the redshifts we are considering. Since we are also resolving Jeans scale in our 
simulation, we can directly compare our theoretical results with this data without taking 
into account the scale dependence of t c q in our analysis. 

The mean opacity (r) and the effective opacity are computed by averaging over all the 
realizations of r for a given redshift bin. 
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i — r 




10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 



k(h ~ 1 Mpc) -> 

Fig. 1. — The matter power spectrum for magnetic case, with added exponential cutoff and 
then smoothed around magnetic Jeans length kj , is displayed for various values of magnetic 
field strength Bq. Spectral index n is -2.95 for each case. Along with that the red curve is 
matter power spectrum for pure ACDM non-magnetic case. 
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5. Results 



In Figure 1 we show the matter power spectra at the present epoch for magnetic case, 
along with the pure ACDM (non-magnetic) matter power spectrum, which has been used 
in our calculations, here an exponential cut-off around k magnetic Jeans scale is assumed. 
This figure shows that the magnetic field induced matter power spectra can dominate over 
the pure ACDM case at small scales (fc > 1 h/Mpc). The effect of this excess has already 



been studied in the context o 



signals ( jSethi fc Subramanian 



early 



(12005 



struct u re formation, rei o nizati on, and weak-lensing 



20091 ). IPandev & Sethi! ( 12012^ As an extension to 



that body of work we explore the effect of this excess on Lya effective opacity in this paper. 

In Figure 2 we show the variation of Lyct opacity (r) with redshift for various values 
of magnetic field strengths. The red dots with y-errorbars are the observed values of Lya 



opacity t c s (IFaucher-Giguere et al. 



(120081 )). It should be pointed out that the inclusion of 
peculiar velocities in the computation of r (Eq. (17)) makes a negligible difference to the 
value of either average or effective opacity. 

Figure 1 corresponds to the case when matter perturbations induced by primordial 
magnetic fields and the inflationary matter perturbations are not correlated. The average 
opacity (r) is not an observable quantity.The aim of Figure 2 is to demonstrate that the 
inclusion of PMF matter perturbations enhances the average opacity of the IGM. 

In Figure 3 we show the variation of t c s with redshift for various values of magnetic 
field strength along with the observed evolution of t c q. This plot is for the case when 
matter perturbations induced by primordial magnetic field and the inflationary matter 
perturbations are not correlated. Comparing this figure with Figure 2 we see that the slope 
of redshift evolution of r c s is far smaller than for average opacity. This difference is owing 
to the fact that for HI column densities Nm > 10 14 cm -2 , the optical depth exceeds one. 
For column densities larger than this saturation value, the optical depth increases only as 
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logarithm of the column density and therefore these clouds get a smaller weight in the the 
computation of r e g. As the average opacity of the IGM increases sharply with increasing 
redshift (Eq. (19)), this effect is more enhanced at higher redshifts. 

A comparison between Figures 2 and 3 shows that an increase in (r) doesn't necessarily 
lead to an increase in r e g. In Figure 3 it is seen that r c ff is greater than the usual ACDM 
case for z < 3 but falls below the predictions of this model for larger redshifts. 

We can understand this behaviour by the following set of arguments. The change 
in the effective optical depth dT c s oc exp(— Ti)dri, where Tj refers to optical depths of 
individual clouds. On the other hand, d(r) oc X]i^ r «- As seen m Figure 2, the inclusion 
of PMF density perturbations increase (r) or J2i^ T i > 0' but £\ exp(— r,j)dTi could be 
negative if dri is negative wherever T{ is smallest. To elaborate this point, In Figure 6 we 
have plotted the distribution of optical depths r^s for the 1.5 nG case [z = 4), against the 
dri = n (2 nG) — Ti (0 nG); (z — 4). It is clear from this figure that Ti values are small 
when dTi is more negative, or this can make £V exp(— r^dr, negative, and thus it explains 
the decrease of r e fj even when there is increase in (r) with increasing magnetic field values. 
It should be pointed out that this behaviour of T c g cannot be mimicked by a change in J, 7 
(Eq. (19)) or by a scaling of the power spectrum by changing the value of as- 

The Figure 4 and Figure 5 are for the same analysis as Figure 2 and Figure 3 
respectively, but for the case when induced matter perturbations and inflationary matter 
perturbations are correlated. In Figure 4 the values of (r) are smaller in comparison to the 
corresponding values in the case of uncorrelated matter perturbations. 

For detaile d comparison with obs e rvatio ns we performed the likelihood analysis for the 



T c ff against the 



Faucher-Giguere et al. 



(120081 ) data as a function of four parameters, J ((1.4 
to 2.0)10~ 12 ), 7 (1.4 to 2.0), B ((0.1 to 2.0) nG), and n (-2.80 to -2.99). To compute the 
posterior probability for magnetic field parameters, we marginalized the likelihood function 
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over the parameters J & 7. Figure 7 shows the results of this analysis for the uncorrelated 
case. The curves from top to bottom are the contours for 5a, 3a and la levels for a range 
of Ax 2 = xf ~ Xmin- We see that in this case for n = -2.90 the allowed values (by 5a) of 
magnetic field are Bq < 0.6-0.7 nG, and for n = -2.95 is -Bo < 1.3 nG. 

n Figure 8, we compar e this result with our previous analysis with the weak-lensing 



data iPandey &: Sethi! (|2012( ) and the present analysis with the correlated case: the lower 
triplet (red green and blue), solid and dashed corresponds to the uncorrelated and the 
correlated cases respectively of the present analysis, whereas the upper triplet (dotted) 
correspond to our previous analysis with the weak-lensing data. It is clearly seen from 
Figure 8 that the constraints arising from the correlated case are not very different from the 
uncorrelated case. Or the Lyman-a clouds do not provide an appropriate physical setting 
for distinguishing between these two cases. From Figure 8 it also follows that our present 
constraints are considerably stronger that our previous analysis with the cosmological 
weak-lensing data. 
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B = 0.0nG 
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Fig. 2.— The evolution of Lya opacity (r) for the magnetic and non magnetic cases, 
uncorrelated S in Q and <5 pm f case. 
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B = 0.0 nG 

B = 0.4nG 

B . 0.6 nG 

B = 1.0 nG 

B = 1.5nG 
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2.5 



3.5 



4.5 



Fig. 3.— The evolution of Lya opacity r cff for the magnetic and non magnetic cases, 
uncorrected <5 in fl and 5 pm f case. 
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6. Discussion 



Primordial tangled magnetic fields leave their signatures on cosmological observables 
for a large range of scales from sub-Mpc to 10 4 Mpc. CMBR temperature and 
polarization anisotropics provide pr obes for the magnetic fields for scales > 10 Mpc (e.g. 



Kahniashvili et al.l 



201(f ). Recently, 



Yamazaki et al 



( 120101 ) computed the allowed region 



in the {Bo,n} plane by comparing the predictions of primordial magnetic field models 
with existing CMBR observations. Constraints on smaller scale come from early formation 



of structures induced by PMF. The observables that impact these scales in c 



reionization, HI signal from the epoch of reionization ( iSethi &; Subramanianl (12009 



Schleicher fc Miniatil (j201lh ). cosmological weak gravitation lensing ( iPandey fc Sethi 



ude early 



2005|) 



( 20121 )) . etc. Other constraints on large scale cosmological m agnetic fields arise 



rotation measure of high redshift polarized radio sources (e.g. 



Kolatt 



1998 



arise 


rom 


Sethi 


2003 



Blasi et al. 



19991 ); rotation measure (RM) of radio sources will be one of the methods 



employed by radio i nterferomete r s LO FAR and SKA to attempt to detect cosmic magnetic 



fields. In particular, 



Blasi et al. 



( 119991 ) considered the same physical setting (high redshift 
Lyman-a clouds) as in this paper. They computed the RM of Lya density field and 
obtained bounds ~ 10~ 8 G on magnetic fields with coherence length scales of the thermal 
Jeans length. 

In addition to the upper bounds on the magnetic field strength obtained by these 
observables, recent results s uggests that t h ere might be a lower bo u nd of ~ 10~ 15 G on th e 



magnetic fie 



Taylor et al 



stren gth (e.g. 



DolaeJ ( I2010h : iNeronov & Vovkl ( 120101 ) 



Tavecchio et al. 



feoid); 



( 1201 ll )). An other lower bound is ob tained from the study of echo emission 



from the blazar Mrk 501 (iTakahashi et al. I ( 120121 )) which suggests magnetic field strength 
of B > 10~ 20 G coherent over the length scale of ~ 1 kpc. This would suggest that the 
magnetic field strength could lie in the range 10~ 20 < Bq < a few 10~ 9 G. This range is still 
too large for a better determination of the magnetic field strength. 
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observation i — i — i 




1.5 2 2.5 3 3.5 4 4.5 

Z — > 



Fig. 4.— The evolution of Lya opacity (r) for the magnetic and non magnetic cases, 
correlated S in Q and 5 pm f case. 
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Fig. 5.— The evolution of Lya opacity r cff for the magnetic and non magnetic cases, 
correlated S in Q and S pm f case. 
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Fig. 6. — The distribution of r, (1.5 nG) versus dri (= 7j (1.5 nG) — r* (0 nC)) at redshift 
z = 4. 
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Fig. 7. — Allowed (t he shaded 1 ) region in th e (Bp, n) plane, based on the \ 2 analysis of r e ff 



against the data from 



Faucher-Giguere et al. 



(120081 ). The three curves (from top to bottom) 



are contours at the 5<r, 3a and la levels. 



21 



T 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

B (nG) ^ 



Fig. 8. — Contour of A% 2 = X\ — Xmm f° r the present analysis (the lower triplets, solid ones 
correspond to 1, 3 and 5 sigma levels for the uncorrelated case, the dashed ones are the same 
for the correlated case) and the previous analysis (the upper triplet) on (B , n) plane. 
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In Figure 7, we show the constraints from the present study comp are with similar 



constr aints from cosmological gravitational lensing we obtained earlier (jPandey fc Sethi 



( 120121 )). In Com parison to bounds on p rimordial magnetic fields from CMBR anisotropics 



(e.g. Figure 1 of 



Yamazaki et al. 



( 120101 )). for the entire range of spectral indices, we obtain 



stronger limits on B . Othe r constraints from b i spect rum and trispectrum analysis of 



CMBR passive scalar modes 



Trivedi et al. 



(12010 



20121 ) are 2.4 nG and 0.7 nG, they have 



used spectral index value n = -2.8, whereas for n = -2.8 our analysis gives an upper bound 
on B < 0.3 nG (Figure 7). As noted above, these bounds are even better than our previous 
analysis with the weak-lensing data (Figure 8). 

In our present analysis, we consider four parameter: J, 7, Bo, n but no other 
cosmological parameters. We also do not account for errors arising from different 
realizations of the density field. Our current bounds can be further improved by the 
inclusion of such effects. We note that even though the magnetic field signal could be 
degenerate with the overall normalization of the ACDM model as measured by as, the 
curre nt errors on the value of erg (WMAP 7-year data give <7g = 0.801 ± 0.030 |l arson et al. 
20111 )) are too small to sufficiently alter our conclusions. 



In sum: Lyman-a clouds provide a sensitive probe of the matter power spectra at scales 
< 1 Mpc. Primordial magnetic field induced matter perturbations give additional power at 
these scale which can be probed using the redshift evolution of r c g. Our results shows that 
this leads to one of the most stringent bounds on the parameters of primordial magnetic 
fields. These bounds can be further improved by more data on th e evolution o f r P ff a t low 



redshifts and also more precise data at higher redshifts. Recently, 



Becker et al 



( 120121 ) have 



provi ded a measurement of the evo lution r c ff which is in agreement with the data we have 



used (IFaucher-Giguere et al. 



(120081 )) for our analysis but claims better precision. In future, 
similar analyses with such data can give even more stringent constraints on the parameters 
of primordial magnetic fields. 
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